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Direct measurement of 7 using B — > Kmr and B — > KKK decays 
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Using the BaBar measurements of the Dalitz plots for B — > , B° — > K°n + n , E> + — > 

K + tt + tv~, B° ->■ K + K°K~, and 5° -> K°K°K° decays, we cleanly extract the weak phase 7. 
We find four possible solutions: (31±l)° , (77 ± 3)°, (2581a) , and (3151a) - In all cases the error 
includes first-order flavor-SU(3) breaking effects. One solution - (77 ± 3)° - is consistent with the 
standard model; its error is far smaller than that obtained using two-body B decays. 
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One of the main aims of B physics is to test the stan- 
dard model (SM) explanation of CP violation, which 
is that it is due to a complex phase in the Cabibbo- 
Kobayashi-Maskawa (CKM) quark mixing matrix. To 
this end^one measures the three angles of the unitarity 
triangle a, /3 and 7, in many different ways, and looks 
for discrepancies. 

The standard method to directly probe 7 was originally 
suggested by Gronau, London, and Wyler (GLW) 
with variants subsequently proposed by Atwood, Duni- 
etz and Soni and by Giri, Grossman, Soffer and Zu- 
pan Q. It involves the measurement of the rates for the 
two-body decays B ± — > DK^, where D = D°, D° or 
D C p = (D° ± D°)/y/2 (D C p decays into a CP eigen- 
state). Although the two-body method is expected to 
be theoretically clean, it is difficult experimentally, so 
that the present direct measurement, using a combina- 
tion of the techniques of Refs. 043 > nas a iar gc error: 
7 = (66 ± 12)° j|. 

Recently it was shown that, contrary to previous think- 
ing, it is possible to extract CKM phase information from 
charmless three-body B decays |g, 0- There are three 
ingredients. First, the three-body decay amplitudes can 
be expressed in terms of diagrams similar to those of 
two-body B decays These are described in Ref. @. 
(As we consider b — > s transitions, the B + decay am- 
plitude can receive a contribution from the annihilation 
diagram. This is neglected.) Note that, unlike the two- 
body diagrams, the three-body diagrams are momentum 
dependent. 

Second, it is possible to fix the symmetry of the final 
state. This is done using the Dalitz plot of B — > P1P2P3 
(the Pi are pseudoscalar mesons) Q . Denoting by pt the 
momentum of each Pi, one defines the three Mandclstam 
variables Sjj = (pi+pj) 2 . These are not independent, 
but obey S12 + S13 + S23 = m 2 B -\-m\ + m>2 + Now, the 
Dalitz plot is given in terms of two Mandclstam variables, 
say S12 and S13. The key point is that the experimental 
Dalitz-plot analysis allows one to reconstruct the decay 



amplitude Ai(B —> PiP2Ps)(si2, S13). The amplitude for 
a state with a given symmetry is then found by applying 
this symmetry to M(si2, S13). This amplitude is used to 
compute all (momentum-dependent) obscrvablcs for the 
decay. For example, the final state Ksir + ir~ has CP + 
if the 7T + 7r~ pair is symmetrized. The amplitude for this 
state is [M(si 2 ,sia) + M(s 13 , si 2 )]/\/2. 

Third, in Rcf. Ml it was shown that, as is the case in 
two-body decays [3], under flavor SU(3) there are rela- 
tions between the electroweak penguin (EWP) and tree 
diagrams for b — > s transitions. These take the simple 
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in which the a are Wilson coefficients and A 



V* h Vn. 



Note: the EWP-tree relations hold only for the state 
that is fully symmetric under exchanges of the final-state 
particles. However, the amplitude for this state can be 
found as described above using the Dalitz plot: 

M is = -^= [M{s 12 ,s 13 ) + M(s 13 ,s 12 ) + M{s 1 2,s 23 ) 
V6 

+ M(s 23 , S12) +M(s 23 , si 3 ) +M(s 13 , s 23 )] , (2) 

where the subscript "fs" stands for "fully symmetric." 

Based on the above, a method was proposed for ex- 
tracting the weak phase 7 from B — > Kttit and B — > 
KKK decays 0. It works as follows. The follow- 
ing five decays are considered: B° — > K + tt°it~, B° — > 
A" 7r+7r-, B+_ -> K+ir+n-, B° -> K + K°K~ , and 
B° —> K°K°K°. In writing the amplitudes of these five 
processes in terms of diagrams, we note the following. In 
three-body decays, one has to "pop" a quark pair from 
the vacuum. For B —¥ Ktttt decays, the popped quark 
pair is uu or dd (under isospin, these diagrams are equal), 
while the B — > KKK decays may have a popped ss pair. 
Now, the imposition of the EWP-tree relations assumes 
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flavor SU(3) symmetry. But this also implies that dia- 
grams with a popped ss quark pair are equal to those 
with a popped uu or dd. In other words, under flavor- 
SU(3) symmetry the diagrams in B — > KKK decays are 
the same as those in B — > Kirn decays. 

Note, however, that flavor-SU(3) symmetry is not ex- 
act. It is therefore important to keep track of a possible 
difference between B — > Ktttt and B — > KKK decays. 

The amplitudes for the five decays in terms of dia- 
grams are given in Ref. ■ We define the following four 
effective diagrams: 



a = -PL 
b = T{ + C 2 , 
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d = T[ + C[ 
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The decay amplitudes can now be written in terms of five 
diagrams, a-d and P' uc : 



2A(B [) -^K+tt^-) 



be n - kc , 



V2A(B° -> K°ir + 7r-)i s = -de 1 ' 1 - P' ac e^ -a + nd , 
V2A(B+ -S- A"+7r+7r-) fs = - ce h - P' uc e n - a + Kb , 
V2A(B° -> K+K°K-) [S = a su(3) x 

(-ce i7 - P' uc e l '< -a + Kb), 
A(B° -> K°K°K°) is = a su(3) (P' uc e^ + a) , (4) 

where ctsu(3) measures the amount of flavor-SU(3) break- 
ing. 

In the flavor-SU(3) limit (\a SU (3)\ = 1), M B+ ~> 
K+it+iT-) is = A(B a -> K+K°K-) is , so that the B+ de- 
cay does not furnish any new information. The remain- 
ing four amplitudes depend on ten theoretical parame- 
ters: five magnitudes of diagrams, four relative strong 
phases, and 7. But there are eleven experimental ob- 
servables: the decay rates and direct CP asymmetries for 
the four B° decays, and the indirect CP asymmetries of 



B° -> K 



B° -> K+K°K- and B° -> K G K Q K [ 



With more observables than theoretical parameters, 7 
can be extracted from a fit. Furthermore, if one allows 
for SU(3) breaking (|as[/(3)l 7^ 1): we can add two more 
observables: the decay rate and direct CP asymmetry for 
the B + decay. In this case it is possible to extract 7 even 
with the inclusion of |asj/(3)| as a fit parameter. 

As has been stressed above, the diagrams and observ- 
ables are both momentum dependent. That is, the above 
method for extracting 7 in fact applies to each point in 
the Dalitz plot. However, the value of 7 is independent of 
momentum, so that the method really represents many 
independent measurements of 7. A preliminary analy- 
sis presented in Ref. considered a naive average over 
such independent measurements of 7. In this letter we 
improve upon this and perform a combined likelihood fit 
to extract 7 from multiple Dalitz-plot points. 

The first step in performing a fit is to collect the ob- 
servables. This is done as follows. An isobar model is 



used to analyze the three-body Dalitz plots. Here the de- 
cay amplitude is expressed as the sum of a non-resonant 
and several intermediate resonant contributions: 



M(s 12 , si 3 ) = A/bp Cje l ^F 3 {s 12 , s 13 ) 



(5) 



where the index j runs over all contributions. Each con- 
tribution is expressed in terms of isobar coefficients Cj 
(amplitude) and 9j (phase), and a dynamical wave func- 
tion Fj. A/dp is a normalization constant. The Fj take 
different forms depending on the contribution. The Cj 
and 9j are extracted from a fit to the Dalitz-plot event 
distribution. Such fits were performed by BaBar [l2| for 
each of the five three-body decays of interest. The iso- 
bar coefficients found, together with their assumed wave 
functions (Fj). allow us to reconstruct the amplitude 
for each three-body decay as a function of the relevant 
Mandelstam variables. We have chosen the normaliza- 
tion constant such that the integral of \A4\ 2 over the 
kinematically-allowed Dalitz-plot phase space gives the 
experimental branching fraction (Bsxp)- We then con- 
struct Aif s using Eq. ([2|). This process is repeated with 
the information available for the CP-conjugate process, 
where we construct Aif s . 

The experimental observables are then obtained as fol- 
lows: 

A(si2,si 3 ) = |A4 fs (s 12 ,s 13 )| 2 + pW fs (si 2 ,si 3 )| 2 , 
Y(s 12l s 13 ) = |A4f s (si2,si 3 )| 2 - |A7 fs (si2,si 3 )| 2 , 
Z(si 2l s u ) = Im [A4f 8 (si2,si 3 ) Atf s (si2,si 3 )] . (6) 

In order to obtain the experimental error bars on these 
quantities we allow the input isobar coefficients to vary 
over their lcr allowed ranges. Note that, although the 
errors on the isobar coefficients extracted from a given 
Dalitz plot are in general correlated, such information is 
not always publicly available. In our analysis we have 
considered the errors to be completely uncorrelated, but 
we hope that a future analysis by an experimental col- 
laboration will take such effects into account. 

The effective CP-averaged branching ratio (X), direct 
CP asymmetry (Y), and indirect CP asymmetry (Z) may 
be constructed for every point on any Dalitz plot. How- 
ever, when the final state has a specific flavor, such as 
in the case of B° — > K + ir°Tr~ , the quantity Z has no 
physical meaning and is therefore left out of our anal- 
ysis. In addition, note that, since the amplitudes used 
to construct these observables are fully symmetric under 
the interchange of the three Mandelstam variables, for 
any given point on a Dalitz plot there will be five other 
points where the extracted X, Y and Z take identical 
values, and hence do not provide any new information. 
In order to avoid counting the same information multiple 
times, we therefore divide each Dalitz plot into six zones 
by its three axes of symmetry, and use information only 
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FIG. 1: Kinematic boundaries and symmetry axes of B —¥ 
Knn and B — > KKK Dalitz plots. The symmetry axes divide 
each plot into six zones, five of which are marked 2-6. The 
fifty dots in the region of overlap of the first of six zones from 
all Dalitz plots are used for the 7 measurement. 

from one zone. This is shown in Fig. [TJ where we select 
the dotted zone for our calculations. 

Although it is possible to measure both the direct and 
indirect CP asymmetries in B Q — > KsK$Ks, their mea- 
surement is currently statistics limited. The experimen- 
tal Dalitz-plot analysis done by BaBar makes no distinc- 
tion between the amplitude and its CP conjugate. That 
is, they take A(B° K S K S K S ) = A(B° K S K S K S ). 
This has two consequences. First, Y and Z vanish for 
every point of the Dalitz plot. Second, this requires that 
P' uc be set to zero in Eq. (U]). The removal of an equal 
number of unknown parameters (amplitude and phase of 
P^ c ) and observables does not affect the viability of the 
method described above. 

With the observables in hand, we now perform a max- 
imum likelihood analysis for extracting 7. We select 
points within a sixth of the overlapping kincmatically- 
allowed region of the Dalitz plots separated from one an- 
other by 1 GeV 2 in both Mandclstam variables. We find 
that there are fifty such points. (Note that the 1 GeV 2 
resolution is ad hoc. In a more appropriate analysis, one 
should consider the optimal bin size for the Dalitz plot 
with the lowest available statistics to determine this res- 
olution.) For each of these fifty points we construct the 
— 2AlnL(7) function, where L represents the likelihood, 
which we then minimize over all the hadronic parameters 
for that point. The sum of such functions over all fifty 
points gives us a joint likelihood distribution. The local 
minima of this function are then identified as the most- 
likely values of 7. In order to find the la error bar on 
7 we first shift the likelihood function along the vertical 
axis so that the zero of the function corresponds to a local 
minimum. We then look for the range of 7 that results 
in a unit shift along the vertical axis of the — 2 A In 1,(7) 



vs 7 plot. 

We perform the likelihood maximization fit in three 
different ways and plot our results in Fig. [2] We first 
consider the scenario in which flavor SU(3) is a good 
symmetry. That is, we fix |o;5(7(3)| = 1; our analysis 
involves only the four B° decay channels. The most- 
likely values of 7 obtained in this way are listed under 
Fit 1 in Table[U 

Second, we allow for SU(3) breaking and treat it as 
follows. We compare the Dalitz plots for the two pro- 
cesses B+ K+Ti+TT- and B° K+K a K~ point by 
point. Theoretically, the amplitudes for these processes 
differ only by the parameter a su ^ . The ratio of Vs con- 
structed from the two Dalitz plots then gives us \oesu(3) \ 2 ■ 
(Note that a similar ratio constructed from the Vs has 
an enormous error due to the smallness of Y . We are 
therefore unable to extract any interesting physical in- 
formation from such a ratio.) Averaged over the fifty 
points we find lag^)! = 0.97 ± 0.05. This shows that, 
on average, SU(3) breaking is small. We use |as[/( 3 )| 
found in this way to correct the observables from the 
B —J- KKK Dalitz plots and use the corrected numbers 
in a new maximum-likelihood analysis for finding 7. We 
present the results under Fit 2 in Table [I] 

In the third maximum-likelihood analysis, we consider 
observables from all five Dalitz plots but now include 
l a st/(3)| as an additional unknown hadronic parameter. 
The results from this method are listed under Fit 3 in 
Table D 

TABLE I: Most likely values of 7 (in degrees) extracted from 
Fig. [2] Results are obtained using the three different fitting 
methods as explained in the text. 



Solution 


Fit 1 


Fit 2 


Fit 3 


I 


3113 


31 ±2 


311* 


II 


76l^ 


781^ 


77 ±3 


III 


26lll 


259 ±3 


2581^ 


IV 


314 ±2 


315 ± 2 


3151^ 



The maximum-likelihood analysis indicates that, in 
each of the three methods described above, the data favor 
four distinct discretely-ambiguous values of 7. In Table 
U we present the most-likely values of 7 extracted using 
these three methods. It is evident from the results that 
the inclusion of an SU(3)-breaking parameter |asc/(3)| 
shifts the preferred values of 7 by only a tiny amount. 
This indicates that the leading-order effects of flavor- 
SU(3) breaking are well under control in three-body B 
decays. While one cannot completely remove this source 
of theoretical error from our analysis, the error bars are 
rather small. 

Even though there are four preferred values of 7, in all 
cases the error is small, 2-4°. Although this may seem 
surprising at first sight, it really is not when one remem- 
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FIG. 2: Results of maximum-likelihood fits. The solid (black) curve represents the fit assuming flavor-SU(3) symmetry. The 
short dashes (red) represent the fit where fiavor-SU(3) breaking is fixed by a point-by-point comparison of Dalitz plots for 
B + — > K + tt + iy~ and B° — s> K + K° K~ . The long dashes (blue) represent the fit with inputs from five Dalitz plots and an extra 
hadronic fit parameter |as^(3)|. 



bers that there are, in fact, fifty independent measure- 
ments of 7. Roughly speaking, if each measurement has 
an error of ±20°, which is somewhat larger than other 
methods, then when we take a naive average, we divide 
the error by VoO, giving a final error of ~ 3°. The one 
caveat is that all errors considered to this point have 
been entirely statistical - the systematic error has not 
been included (we do not know its value). Hopefully, 
the experimentalists themselves will redo this analysis, 
including all errors. 

There are other sources of error that have not been 
included in our method. First, as mentioned above, we 
have assumed that the errors on the isobar coefficients 
are uncorrelated. In practice such errors are correlated, 
which will alter the error bars on X, Y and Z . Second, we 
have only taken leading-order flavor-SU(3) breaking into 
account. Higher-order flavor-SU(3) breaking may arise 
due to the nonzero mass difference between pions and 
kaons, and between intermediate resonances. This said, 
the error due to leading-order SU(3) breaking is evidently 
small. It is unlikely that the error due to higher-order 
SU(3) breaking is larger. 

Finally, we reiterate that we find four most-likely val- 
ues of 7 in B — > Ktttt and B —¥ KKK decays. Three 
of these - (31±3)°, (258±3)°, and (315j^)° - are in dis- 
agreement with the SM (is this a "Kirir-KKK puzzle"?). 
However, one solution - (77 ±3)° - is consistent with the 
SM. Furthermore, its error (which includes SU(3) break- 
ing, but not other systematic effects) is far smaller than 
that found with the GLW method, which uses two-body 
B decays. The fundamental reason is that three-body de- 
cays give many independent measurements of 7. These 
can be combined, reducing the error on 7. 
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